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New physics scenarios beyond the Standard Model predict the existence of milli-charged particles. 
So far, only spin-1/2 and spin-0 milli-charged particles have been considered in literature, leaving 
out the interesting case of spin-1. We propose a minimal unitary and renormalizable model of 
massive milli-charged vector particles. Unitarity requires that these particles are gauge bosons of 
a non-abelian spontaneously broken gauge symmetry. The minimal scenario then consists of an 
extended Standard Model gauge group SU{2)l x U{1)y x SU{2)d together with a SU{2)d dark 
Higgs boson responsible for the symmetry breaking in the dark sector. By imposing that the dark 
Higgs multiplet has a non-vanishing milli-hypercharge, stable milli-charged spin-1 helds arise thereby 
providing a potential dark matter candidate. We analyse the phenomenological constraints on this 
scenario and discuss their implications. 


I. INTRODUCTION 

Although we observe electric charge quantization in 
Nature, this property is not a requirement of the Stan¬ 
dard Model (SM) [T]. New physics beyond this frame¬ 
work could enforce the charge quantization, grand uni¬ 
fied theories [5] for instance, and the existence of mag¬ 
netic monopoles, if experimentally corroborated, would 
demand it |3]. 

Other theories beyond the SM, however, predict the 
existence of particles with a small, either non-quantized 
or effective, electric charge ee. Given that e ^ 1 and e 
is the electron charge, these particles are usually referred 
to as milli-charged particles (MGP) [IHlj. MCP arise 
naturally in theories where two or more unbroken C/(l) 
gauge symmetries, coupled to different matter sectors, 
possess non-diagonal kinetic terms. Even in the absence 
of a tree-level kinetic mixing, a non-diagonal kinetic term 
is inevitably induced by radiative corrections [5]. 

MCP are stable because of charge conservation and 
thus are a natural candidate for dark matter. In this 
regard, their direct coupling to the SM photons also 
provides suitable production mechanisms as the thermal 
freeze-out mi or freeze-in m- The possibility of vec¬ 
tor MCP has not yet been explored in literature, where 
the dedicated studies focused on the contexts of spin-0 
and 1/2 [TT1[T2]. However, the case of spin-1 MCP, V/, 
presents and intriguing feature: as a result of the inter¬ 
play between gauge interactions and unitarity, the total 
cross section of 77 —)■ VV tends to a constant in the high 
energy limit, whereas the same quantity decreases as 
in the cases of spin-0 and spin-1/2 MCP. Such distin¬ 
guishing characteristic of spin-1 interactions is manifest 
in the SM, where the tree-level total cross section for 
77 —> W^W~ approaches a constant of about 80 pb at 
high energies [I3HI1], while radiative corrections are typ¬ 
ically of order 10% [16] . 

Furthermore, while the interactions of scalars and 
fermions with the photon can always be induced via ki¬ 


netic mixing of different U{1) gauge sectors, the spin-1 
case requires the extension of the SM gauge group to a 
larger non-abelian gauge group. Motivated by these ap¬ 
pealing characteristics, we present here a minimal model 
of spin-1 milli-charged particles, which preserves both 
unitary and renormalizability, and investigate its phe¬ 
nomenological implications. 

Adopting a minimalistic approach, we consider the 
SU{2)l X U{\)y X SU{2)d group, where SU{2)d is the 
gauge group of the dark sector that operates on the com¬ 
ponents Vfi of a dark gauge multiplet. In order to provide 
the latter with a mass, we also introduce a SU{2)d Higgs 
multiplet with a non-zero hypercharge. After symmetry 
breaking occurred in both the dark and visible sectors, 
two of the vector fields acquire masses and milli-charge 
couplings e to the ordinary photon. The remaining vec¬ 
tor field also acquires a mass, however it does not couple 
to the photon. The proposed model then resembles the 
Stiickelberg Z' extension of the SM HMS], with the im¬ 
portant difference that a non-abelian gauge group is at 
the basis of our construction. 

In the absence of SU (2) d matter fields other than the 
dark Higgs multiplet, the conservation of electric charge 
guarantees that our massive vector MCP are stable and, 
therefore, a potential dark matter candidate. In this 
regard, the model of spin-1 MCP we present naturally 
yields a rich and interacting dark sector without intro¬ 
ducing a dark or hidden photon [201125] . Whereas the 
phenomenological bounds on the latter clearly do not 
affect our construction, in the following we will inves¬ 
tigate the viability of the proposed scenario by analysing 
the relevant collider, astrophysical and cosmological con¬ 
straints. 


The paper is organised as follows. In the next Sec¬ 
tion we present the most minimal renormalizable model 
for charged spin-1 fields, in Section HI we will analyse 
the phenomenology of this scenario, ranging from astro¬ 
physics and cosmology to low-energy experiments. Our 
conclusions are reported in Section EYI 
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II. MILLI-CHARGED VECTOR BOSONS AND 
THE STANDARD MODEL 

A. Lagrangian for charged spin-1 fields 

Charged and massive spin-1 particles can be described 
by complex vector fields V^, /r S {0,1,2,3}. Restricting 
ourself to operators of dimension 4 or lower, the most gen¬ 
eral, U{1)em invariant and parity conserving Lagrangian 
for interacting charged spin -1 fields is [ 2 ^ 

Cem+v = - + ml\V\^ 

- iQv{g - l)F^''VX 

-^(Ai(|Rn"-A 2 |R 2 | 2 ), ( 1 ) 

where Qy is the electric charge of the particles associated 
to V^, g their gyromagnetic ratio and my their mass. 
With F^i, = we denot^the usual field strength 

of electromagnetism, while is the field tensor of the 
charged vector fields after imposing the minimal sub¬ 
stitution: = {d + iQvA)[^V^]. The parameters Ai 

and A 2 govern the self interactions of the fields 

We remark that for fields of spin higher than 1/2, the 
interactions with the photons are not uniquely deter¬ 
mined by the minimal substitution, because of the non¬ 
commutativity of the covariant derivatives. Requiring 
the unitarity of the theory, however, further reduces the 
set of free parameters to the mass my and charge Qv- 
To be precise, the unitarity of 77 —>■ VV scattering sets 
the gyromagnetic ratio to g = 2, while the unitarity of 
VV —t VV necessarily yields Ai = A 2 = 1 [21] . 

In spite of that, the presence of an explicit mass term 
still leads to unitarity violations, although in softer way, 
and spoils the renormalizability of the model. These 
problems are solved by noticing that, if the mass term 
was to be excluded from the Lagrangian in Eq. Q, the 
choice of parameters Ai = A 2 = 1 and g = 2 would en¬ 
large the gauge symmetry of the latter from U{1)em to 
SU{2). By generating the mass terms via spontaneous 
breaking of this enlarged symmetry, SU{2) —>• U{1)em 
for instance, the unitarity and renormalizability of the 
theory then follow [? ]. Notice that the adopted symme¬ 
try breaking pattern does not affect the essence of this 
construction, which relies on the well known result that 
renormalizable and unitary interactions between vector 
fields are gauge interactions [2M25] . 

B. A minimal SM extension with milli-charged 
vector bosons 

In order to accommodate massive milli-charged vector 
fields within the framework of the SM, we extend the 


^ The bracket notation indicates here := — Ai/B^. 
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TABLE 1. The quantum numbers of vacuum in the considered 
SM extension. 


gauge group of the latter to SU (3) x SU (2)l x SU{2)e>^ 
C7(l)y. The particle content of the theory comprises 
three new vector fields, I/^, that transform under the ad¬ 
joint representation of the dark interaction and a scalar 
multiplet of SU{2)d, (f^D, which is responsible for the 
symmetry breaking in the dark sector and possesses a 
tiny coupling to the SM t/(l)y, i.e. a milli-hypercharge. 
As the extended symmetry is broken by the Higgs mech¬ 
anisms in the visible and dark sector to U{1)em, the 
milli-hypercharge of the dark Higgs (j)D induces opposite 
electric milli-charges in two of the dark gauge bosons. 

The SU (2 )D Yang-Mills term for the dark gauge sector 
is 

Fsui2)n = -\v3^.vr - _ ^g^vnv^ 

( 2 ) 

where (?_d is the dark gauge coupling and V^^v = d[^Vi,y] 
and V/y = {d + igDVz)[^V^] are respectively the field 
tensors of the vector bosons V 3 and V = (Vi -|- iVi)l^f2. 
Although the presented expression for the SU{2)e Yang- 
Mills Lagrangian might seem odd, this form makes ex¬ 
plicit that the Lagrangian in Eq. recovers the one in 
Eq. Q if the mass term is omitted, g = 2, Ai = A 2 = I 
and A is replaced by V 3 . 

As previously pointed out, in order to guarantee the 
unitarity of the theory the mass term for the above vec¬ 
tor bosons has to emerge from a Higgs sector. By de¬ 
noting with (f) the SM Higgs doublet and (pD the dark 
Higgs of milli-hypercharge Tqj give the corresponding 
Lagrangian as 

Fh,hD = 

+ \Di,(p\^ - - A0|(/|^ - A^o \(pDt, (3) 

where stands for the covariant derivative 

D^j, = {d- igWWi - ig'BY - igoV^TD^a)^ (4) 

and r, te are the generators of SU{2)l and SU{2)e 
respectively. The extended gauge group is spontaneously 
broken once 4> and (pD acquire vacuum expectation values 
V and U£), respectively. The quantum numbers of the 
vacuum state of the theory are listed in Table |Tj 

C. Milli-charge from mass mixing 

The electric charge of a particle species can be de¬ 
duced from the following two basic facts: first, a residual 
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U{1)em, or equivalently a massless photon, is possible 
only if the vacuum is electromagnetically neutral, and 
second, the electric charge Q of any particle species can 
be expressed as a linear combination of the hypercharge 
Y, the weak isospin J 3 and the dark isospin J 3 D. Given 
the Higgs charge and isospin assignments listed in Ta¬ 
ble |Tj setting the electron charge to —e uniquely fixes the 
charges of the remaining particle content according to 
the generalised Gell-Mann-Nishijima formula 

Q = e{h + el3D+Y), (5) 

where e = —Yq/I^e q. It is then evident that every parti¬ 
cle with a non-vanishing dark isospin carries also an elec¬ 
tric charge. In particular, the three dark gauge bosons 
with dark hypercharges I 3 E = 0 , ±I result after the sym¬ 
metry breaking into a neutral Z' boson and two vector 
bosons V, , with an electric charge ee. 

Let us remark on the structure of the adopted scheme. 
The most common alternative for generating MGP is via 
a kinetic mixing of two or more 1/(1) gauge fields. Apply¬ 
ing such a scheme to the spin -1 case, however, necessarily 
requires that the SM gauge group be extended at least to 
the SU{2)e X U'{ 1 ) grouiQ The proposed model of vec¬ 
tor MCP arising from the mass mixing is, therefore, the 
most minimal choice in terms of gauge and scalar fields, 
containing only renormalizable interactions. 


D. The neutral gauge bosons 

Spontaneous symmetry breaking generates the follow¬ 
ing mass terms for the gauge bosons 

Cm = m^V+V. + m^IP+Ip- -k (6) 

where mw and mv are the masses of the SM and dark 
sector charged gauge bosons respectively. The masses of 
the neutral gauge bosons are instead contained in the ma¬ 
trix M^, which in the basis = (Vs,/?, W 3 ) is written 
as 

( m|, — 0 \ 

-ml^cwsw ■ (7) 

0 -ITT-Zo^WSw / 

Here we introduced the small dimensionless expansion 
parameter 


C = f^9'/9D, ( 8 ) 

to quantify the deviations from the decoupling limit 
e = ^ = 0 in which the SM is recovered. Notice that the 


^ Notice that, up to fields redefinition, imposing a small tree-level 
kinetic mixing among two U(1) X U{!)' gauge fields is equivalent 
to the direct assignment of an effective milli-charge coupling 


mass matrix in Eq. Q matches the one of the Stiickelberg 
extension of SM, while the expansion parameter ^ coin¬ 
cides with the mass ratio Mi/M 2 commonly used in the 
framework of Stiickelberg Z' models [15]. In this study 
we assume e ^ 1 , as suggested by the phenomeno¬ 

logical constraints that the searches of MCP and Z' im¬ 
pose [H]. 

The mass parameters in Eq.s @ and ([7| 

mw = y, rnzo = mw/cw, (9a) 

mv = VD9D\lD3,o\r, mz'^ = my/r, (9b) 

correspond to the tree-level masses in the decoupling 
limit. Here sw, Cfv = 9 /\/ 9 ^ Y g'^ denote the sine and 
the cosine of the SM weak mixing angle 9\y respectively. 
The parameter 



( 10 ) 


quantifies the effect of the dark Higgs vacuum state and 
depends on the representation adopted for the latter. In 
the rest of the paper we consider a SU{2)e dark Higgs 
doublet, corresponding to the most minimal choice, for 
which r = 1 and therefore mz^ = my- The tree level 
neutral gauge boson masses ar e the roots of the polyno- 
mial det(s —M^), given in Eq. (A2) of Appendix [A] with 
the zero mass eigenstate corresponding to the photon. 
The remaining non-vanishing eigenvalues 


rriy = m 


Zo 


Y 




9 9 


o(e"), 


Yeml 


1 - 




2 9 

— ml/ 


( 11 ) 


+ 0(^4). ( 12 ) 


give the Z and Z' boson masses. The simplest approach 
to identify the corresponding eigenvector fields in a per¬ 
turbative fashion is to start from the basis 



0 

Cw 

-Sw 




(13) 


formed by the neutral gauge boson mass eigenstates in 
the decoupling limit. The mass eigenstates of the mass 
matrix in Eq. Q are then given through the rotation 



( C1C2 -S1C2 -S2\ / ^o\ 

Si Cl 0 Ao , 

C1S2 —S1S2 C2 / \^o/ 


(14) 


where Ci and Si stand for the cosines and sines of the 
mixing angles. In the small mixing approximation, ^ <gC 
1 , we have 


61 « -^cw, 


^ m^^sw 


m%, — m 


Zo 


(15) 
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Notice that the photon field comprises here an additional 
contribution from Zq but not from Zq, the Z boson in 
the decoupling limit. As a consequence the elementary 
charge in (§ is modified as 

e = ( 5 '”^!+ +5"^) ' «eo-e|]^V2, (16) 

where Cq = e^=o = g'cw- We remark that the leading 
contribution to the charge and masses is of the second 
order in whereas the mixing in Eq. (14) contains also 


first order corrections in this expansion parameter. 

Given the above relations, we choose to parametrise 
our scenario with the following physical quantities: mw, 
mvI a, e, These will serve as input parameters in the 
computations of the observables connected to the phe¬ 
nomenology of vector MCP discussed below, where the 
r parameter in Eq. (10) is set to r = 1 as previously 


explained. Throughout the following we will also ap¬ 
proximate most of the quantities with the lowest order of 
their ^ expansion. The subscript that denotes quantities 
in the decoupling limit can therefore be omitted, with 
the understanding that the committed error is of higher 
order in 

To summarise our results so far, we identified a new 
mechanism that yields massive vector MCP from a non- 
Abelian extension of the SM. Because of charge conser¬ 
vation, our MCP are stable and therefore provide a suit¬ 
able dark matter candidate. Our construction is based 
on the presence of a dark Higgs field, required by the 
unitarity and renormalizability of the model, that trans¬ 
forms non-trivially under the gauged symmetry of the 
dark sector. Provided the dark Higgs has a non vanish¬ 
ing milli-hypercharge, the symmetry breaking dynamics 
then result in a dark sector composed of a neutral vector 
boson, Z', and two further spin-1 particles with oppo¬ 
site milli-charge, . The same dynamics also predict 
modifications to the Electroweak sector of the SM, which 
however are negligible at the current level of precision 
achieved by the dedicated experiments. In order to show 
that the proposed framework is consistent with dark mat¬ 
ter, we now turn our attention to the phenomenology of 
MCP. 


III. PHENOMENOLOGY 

Although MCP can easily avoid the bound from direct 
detection and collider experiments owing to their small 
coupling with the photons, in the past decades dedicated 
experiments as well as astrophysical and cosmological ob¬ 
servations pnUHn] have severely constrained the available 
parameter space. 

More in detail, for light MCP characterised by a mass 
scale rrie well below the electron one, production pro¬ 
cesses like the pair production in an external static mag¬ 
netic fields [ini |4T] exhibit non-perturbative effects that 
lead to large enhancements in the production rates. The 
presence of MCP can then be detected by investigating 


the birefringence and dichroism of polarised laser beams 
that propagate in a strongly magnetised vacuum [38]. 
To date these experiments cast the most severe labo¬ 
ratory constraint on light MCP: e < 0(10“® — 10“"^). 
Strong bounds are also brought by the invisible decay of 
orthopositronium |42] and from the Lamb-shift measure¬ 
ments, which limit the MCP contribution that adds to 
the QED expectations [13]. On the explored mass range, 
10“^ eV < TOe < 10® eV, these experiments impose e < 
0(10“^). From the observational point of view, the pro¬ 
duction of sub-eV MCP in photon-photon collisions yield 
distortions in the energy spectrum of the Cosmic Mi¬ 
crowave Background radiation (CMB) [33133]. The non 
detection of such distortions then implies a bounds of the 
order of e ^ 10“^ on the milli-charge. Stronger but model 
dependent bounds on this parameter are also brought by 
stellar evolution, which constrains e < 0(10“^"*^) on the 
range < 0(10 KeV), and by Big Bang Nucleosynthe¬ 
sis, for which e < 0(10“®) if rUe < 0(10MeV) [44] . 

Heavier MCP candidates have been investigated in ac¬ 
celerator experiments, which resulted in e < 10“"*, 10“® 
respectively for « C>(1) MeV and We « 0(100) MeV 
at the 95% confidence level [30]. A current proposal for 
a new experiment at the LHC, |45j . aims to constrain 
MCP with masses 0.1 < < 100 GeV and couplings 

10“® ^ e < 10“^. On the cosmology side, strong bounds 
for rUg < 10 GeV have been derived by analysing the 
impact of MCP on the angular power spectrum of CMB. 
MCP that are kinetically coupled to the baryon and elec¬ 
trons at the recombination era take part in the acoustic 
oscillations of the baryon-photon plasma [33113], result¬ 
ing in contribution to the power spectrum degenerate 
with that of baryons. Independent measurements of the 
abundance of the latter, from BBN for instance, then 
cast a severe upper bound on the abundance of MCP at 
recombination. In particular, for e > 10“®, MCP can 
only constitute a subdominant component of dark mat¬ 
ter. Further studies of the recombination epoch dynam¬ 
ics yield the upper bound e < I0“® for MCP with mass 
me « 1 GeV, which softens to e < 10“^ for heavier par¬ 
ticles with masses of order 10 TeV [13. Whereas such 
particles could be observed at dark matter direct detec¬ 
tion experiments [48] . the corresponding bounds are not 
directly applicable because our galactic magnetic field 
tends to expel the MCP from the galactic disk [13 09] . 

In the following we analyse the most crucial of these 
bounds within the proposed framework of vector MCP. 


A. Z' phenomenology 


Our scenario predicts the existence of a Z' boson, 
which is not stable in our framework. The width of the 
Z' boson is given by 


Sla^^mv 

48c^ 


2 X 10 niv, 


(17) 





5 


Process 

Amplitude order 

w+w-,f+ f - — v+y- 

€ 

■yZ',ZZ' ^ v+y- 


yj,jZ,ZZ ±^V+V- 


z' z' ^ u+ y- 

^/e 

z' ^ f+ f-,w+w- 

€ 


TABLE II. Processes that regulate the abundance of and 
Z' in the early Universe and the relative amplitudes order, 
with respect to their leading dependence by the parameter ^ 
and milli-charge e dehned in section II-D. 


where we used mz' ~ rriy that holds barring correction 
of order see Eq. (12) and the discussion following 
Eq. (9a). 

Our Z' shares many properties with the massive vec¬ 
torial particle appearing in many U(l)' extensions of the 
SM. The negative results of experimental Z' searches can 
then be used to constrain the parameter space of our 
model. 

Notice that, since the mass matrix of Eq. 0 coin¬ 
cides with that of the Z' Stiickelberg extension of the 
SM, both the theories yield the same Z — Z' mixing phe¬ 
nomenology. As a consequence, the experimental bounds 
that constrain the Stiickelberg Z' model can be directly 
applied to ours. In particular, by requiring the compat¬ 
ibility of our framework with the electroweak precision 
tests we obtaii0^ < 0.06 [50] . 

On the othernand, the searches for a Stiickelberg Z' 
require that Mz' > 890(540)GeV for f = 0.06(0.04) [50] . 
From these results, it follows that less restrictive lower 
bounds on our Z' mass can be achieved by decreasing 
the value of the parameter 


B. Dark matter relic abundance 


ii) > 2mz', barring (pjy as a dark matter candidate 

because of the SU{2)d interactions. 

hi) The 4>d-4>SM mixing is negligible. 

iv) e ^ , barring the Z' as a possible dark matter 

candidate and effectively decoupling the dynamics 
of Z' and (see Table [h]) . 

Under the mentioned conditions emerges as the 
only viable dark matter candidate. The stability of these 
particles is insured by the conservation of the electric 
charge, in absence of lighter, coupled, fermion fields of 
charge e. Barring the production via scalar mixing, the 
abundance of is modified only by the s-channel pair- 
production/annihilation from and into the SM charged 
particles and photons. A direct calculation, reported in 
Appendix reveals that the cross section for the 
annihilation into fermions is always larger than the one 
for annihilation into by about two order of magni¬ 
tude. We will then disregard the impact of the latter in 
our computations. The abundance is regulated by 
the Boltzmann equation 011 

nv + 3Hnv = - ^ {crvM 0 i) {jiy - Uy^^q) , (18) 


where a dot stands for the differentiation with respect 
to the coordinate time, {uvm^i) is the thermally aver¬ 
aged interaction rate and ny^Eq is the equilibrium den¬ 
sity of the milli-charged bosons. In the absence of a 
charge asymmetry the total number density of vector 
milli-charged particles is ny := ny+ + ny- = 2ny±. 
Assuming the conservation of entropy, as well as that the 
relevant dynamics takes place in the radiation dominated 
regime, Eq. (181 can be recast as 


da; 


-A(y 2 



(19) 


We investigate now the conditions for our vector MCP 
to populate the dark sector of the Universe. In order to 
identify a tentative dark matter candidate amongst the 
available particles, we first consider the interactions that 
regulate the abundances of V^, Z’, and (/)_d in the early 
Universe. The processes shown in Table |ll] link these 
particle species in a way that the evolution of the cor¬ 
responding abundances can only be tracked by solving a 
set of coupled Boltzmann equation. Given the purpose 
of the present paper, we choose to defer such a detailed 
analysis to a dedicated work |51j and propose here a sim¬ 
pler scenario that nevertheless demonstrates the viability 
of the model. The simplifying assumptions that we con¬ 
sider are 

i) e <C 1, ^ <C 1, hence we disregard the processes with 
amplitudes of higher order in e and 


where Y := ny j s is the comoving density normalised to 
the entropy density s, Yeq := Ueqls, x := m/T and we 
defined 


1/2 / TI" , , 1 


^myM^x^{avM0\) X ^ 


A := g.. 
with the Mpller velocity 

V(Pi ■ P 2 ) - mjml 


VM0I — ' 


E1E2 


( 20 ) 


( 21 ) 


and the labels ‘1’ and ‘2’ are referred to the initial state 
particles. In the above equation g, quantifies the effective 
number of relativistic degrees of freedom and Mpi is the 
Planck mass. The present abundance of spin-1 MCP is 
given by 


riyh^ 


2.8 X 10 


11 my 

TeV 


Mo. 


( 22 ) 


® The ratio Mi/M 2 of 1501 coincides with our parameter 


where Yq, the comoving density of MCP at the present 
time, is obtained by integrating Eq. (19) until Xg = 
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mvIT q, being Tq the CMB temperature of today. 

The parameters of the proposed model are then clearly 
constrained by the requirement that the abundance of 
MCP does not exceed the measured dark matter one: 

< VloMh^ = 0.1199(27) [52]. In order to calculate 
Q.V, we consider two complementary production mecha¬ 
nisms: the thermal freeze-out and freeze-in. 


1. Production via freeze-out 

In our scenario, the relevant process for the freeze-out 
mechanism is the s-channel annihilation of V to SM 
particles. As remarked before, the t-channel processes are 
indeed suppressed by a higher order of the milli-charge. 
The corresponding cross-sections at the threshold s « 
Amy are 


crv+v--^SM — 


16C^ TOy 


387ra^ i 

9 rriy ^ 


(23) 


where (3v = a/I — Arriy/s is the speed of in the cen¬ 
tre of mass frame. The second term corresponds to the 
leading order contribution from V~^V~ —>■ 77 that we re¬ 
port for completeness. Notice that although suppressed 
by a higher power of the milli-charge, this term presents 
a 1 /jfy enhancement at low energies with respect to the 
first one. Hence, although the corresponding contribu¬ 
tion is certainly negligible during the freeze out dynam¬ 
ics, it could give rise to observable effects at later eras 
which we plan to investigate when dealing with the full 
model. 


As the temperature drops below the mass of the par¬ 
ticle, the inverse annihilation rate becomes smaller than 
the expansion rate of the Universe and Eq. (191 reduces 
to Y' = —XY^. Integrating this equation yields a final 
abundance which is inversely proportional to the cross- 
section: 


1 _ 1 

^0 - 



(24) 


Here Xf is the freeze out temperature, determined 
through the condition {YXj^y{xf) = (5(i5-|-l)A(a:/) where 
i5 is a numerical constant of order 0{1). 

The abundance at the freeze out, Yf := Y{xf ), is typ¬ 
ically much larger than the present abundance lb and 
the term Y^^ can then be safely neglected at the desired 
accuracy order. Then, by approximating the thermally 
averaged cross section with the leading order of its low 
temperature expansion 


(o' UM01) 


277ra^ 

16 cfy mlr^ 


(25) 


we obtain the freeze-out temperatur^ 

corresponding to a relic abundance of 

^Mcph^ « 2.4 X lO^x} il^) ■ 

Imposing the upper bound O.MCph'^ < 0.12 results in 
the following rough upper bound on the MCP mass in 
the region e = 0 ( 10 “^ — 10 “®): 

mv<ex 100 GeV. (28) 

Given the constraints mentioned at the beginning of the 
Section, the freeze-out mechanism does not allow for vec¬ 
tor MGP dark matter in this simplified scheme: dark 
matter candidates in the allowed mass range yield an 
overabundance that cannot be depleted in this simplified 
scheme. The necessary higher annihilation cross section 
could however be achieved by considering other annihi¬ 
lation channels, that could dominate over the proposed 
milli-charge channel, or by introducing resonances in the 
latter. The minimal model we propose already contains 
suitable candidates to implement both the mentioned en¬ 
hancement mechanisms and, in a future paper m , we 
intend to study the impact of the Higgs portal and the 
dark Higgs boson as well as the possible role of the Z' 
boson in determining the dark matter relic abundance. 

For instance, allowing for e ~ ^ and therefore gu cx g' 
yields a scenario in which the abundances of and Z' 
are in equilibrium with each other in the early Universe. 
The decays and inverse decays of the Z' boson would then 
provide additional thermal contact between the visible 
and dark sectors and control, along with the reaction 
Z' Z' ^ U+ V ~, the DM freeze-out production. 


2. Production via freeze-in 


The freeze-out scenario is not successful because of the 
too large annihilation cross-section that it requires. It is 
then plausible that the desired spin-1 MCP abundance 
could instead arise via the freeze-in mechanism |10] . In 
this scheme, the particles being produced never achieve 
the equilibrium density, so that Y ^ Yeq and the Boltz¬ 
mann equation in Eq. (191 reduces tcj^ 


Y' = XY'i 


(29) 


^ In our estimates we take the value of the fine structure constant 
as given at the electroweak scale. 

5 Given the relation = ^SM->v+V -'''sm e, 

that holds between the equilibrium densities, we employ the an¬ 
nihilation cross section in Eq. ( | 30 | l to describe the inverse anni¬ 
hilation rate. 
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and the MCP production takes place at energies much 
higher than the corresponding mass scale. The MCP 
equilibrium density is then approximately constant, 
Ygq « 0.28gv/g*, with gy = 6 denoting the degrees of 
freedom of V. The inverse annihilation that drives the 
MCP production is effectively Bose-enhanced, resulting 
in a contribution at most of order 0 ( 1 ) that we however 
neglect in the present analysis. 

For T > my, the relevant cross-section is approxi¬ 
mated as 


o’y+y-->SM — 


StTO^ 2 
16c^s^ 


327ra^ 4 

Only 


(30) 


where, on top of the s-channel fermion annihilation con¬ 
tribution proportional to , we considered the t-channel 
processes 77 —The possible relevance of the 
latter stems from its constant behaviour in the high en¬ 
ergy limit: even if this process is of higher order in e, the 
corresponding cross section does not present the typical 
1 /s suppression at high energies, therefore this channel 
dominates the production as long as the reheating 
temperature is considerably higher than my. 

The corresponding thermal average is then given by 


TTa 


o'VM0i} « Tpe + 


32 


128c4 


2 2 
e X 


(31) 


and by supposing that the MCP have a vanishing initial 
abundance and are produced only via the inverse annihi¬ 
lations, we obtain a relic abundance of 


flh^ « 2.1 X 10^® (e^Xs + 100, (32) 


where Tj. := myjxr is the reheating temperature after 
Inflation and Tg := myjxg is the effective temperature 
at which the production of vector MCP stops. The above 
result holds under the assumption that ^ Tg and 
numerical solutions of Eq. (19) give the value Xg « 1.36. 
The requirement < 0.12 then leads to the following 
minimal upper bound on the milli-charge: 


< min(2 X 10"^°, 5 x 10 


-6 




(33) 


The Xt independent bound, the first term in Eq. (30), 


holds for MCP of any spin within an accuracy of order 
0(1). The second bound, instead, is due to the char¬ 
acteristic behaviour of the spin- 1 —J-spin-l cross section. 
Given that the constraints on the Stiickelberg Z' set 
wy > 0.5 TeV, it follows that the reheating tempera¬ 
ture required to create a DM abundance of vector-MCP 
via this channel is of the order of Planck mass and there¬ 
fore this channel is currently disfavoured by the CMB 
observations. 


C. Astrophysical plasma effects 

The proposed milli-charged dark matter forms a pair 
plasma which can affect the dynamics of galaxy cluster 


collisions by forming shockwaves even if the plasma is 
effectively collision-less [53]. Constraints for dark mat¬ 
ter self-interactions have been derived from observations 
of cluster collisions |S31[SS]. The cleanest example, and 
probably the most constraining, is the Bullet Cluster for 
which it was shown that no more than 30% of the to¬ 
tal DM mass can be stripped from the sub-cluster as it 
passes through the main cluster halo [51] . 

The dynamics of non-relativistic collisions of pair plas¬ 
mas in the non-linear regime, where we expect a fully 
formed shockwave, are still not well understood. On the 
contrary, in the linear regime it is possible to estimate 
the growth rate of plasma instabilities and therefore also 
roughly the time scale of shock shock formation [56] 


10V-0.1yrx(,^) 


10 


-3 


(34) 


1 /2 

where Up = e {At: anjmy) ' is the plasma frequency, n 
denotes the number density of the plasma constituents 
and we assumed a DM density of puM ~ 0.1 GeVcm“®. 
The time scale of a galaxy cluster collision is of the order 
of 0.1 Gyr. If we require that no shockwaves are formed, 
then all dark matter may be milli-charged if 


e < 10 "® 



(35) 


We stress that this inequality is indicative, at best, and 
more accurate methods need to be used for deriving a 
reliable estimate. However, the above result seems to 
indicate that the shock behaviour of the plasma can be 
completely neglected in freeze-in scenarios with heavy 
dark matter my > TeV. If the inequality was violated 
as in the presented freeze-out case, then we would expect 
that astrophysical observations could disprove the vector 
MCP DM scenario. 

Finally we remark that milli-charged dark matter also 
interacts with the intra-galactic visible plasma. This 
might have a non negligible impact on the distribution of 
the X-ray emitting astrophysical plasmas, observed, for 
example, in cluster collisions. A detailed examination of 
these effects is however beyond the purpose of the current 
study. 


D. Acoustic plasma oscillations at the 
recombination 


As previously mentioned, MCP that are tightly cou¬ 
pled to the SM particles at recombination participate in 
the acoustic oscillation of the particle plasma at recom¬ 
bination, yielding a contribution to the CMB spectrum 
that is degenerate with the one brought by baryons. In¬ 
dependent measurements of the abundance of the latter 
then allow MCP to be only a sub-dominant component 
in the observed dark matter relic density. The tight- 
coupling condition, that forbids dark matter MCP, is 










quantified in the following inequality [311 SSI 


> 5 X 10-“- 


mv 




GeV 


- 1/2 


(36) 


where ■= rayirix/{my + Tnrix)- The resulting bound 
on MCP models is reported in Figure where the grey 
area denotes the region of the MCP parameter space in 
which the interactions between MCP protons and elec¬ 
trons are substantial enough to ensure the tight coupling 
of the former to the latter. Given the constraint cast 
by Z' searches, this bound is clearly consistent with the 
candidate we propose. In the above derivation, the in¬ 
teraction between photon and MCP have been neglected 
since the relevant diagrams are of the fourth order in e. 


where we assumed Tyij. « 10~^my and py « 

0.1 GeV cm”^ in order to obtain a conservative estimate. 
As cooling through Compton scattering requires even 
larger timescales 


io 


135 


lompton- 64^3 





.-4 


(39) 

we conclude that DM haloes composed of our spin-1 MCP 
do not collapse to disks for the considered values of the 
parameters. 


IV. SUMMARY 



m\/ (GeV) 


FIG. 1. The exclusion area cut in the MCP parameter space 
by the requirement that MCP do not take part in the acoustic 
oscillations of the particle plasma at the recombination era, 
[37l |46l 113 . In our model, the region my < 500 GeV is 
excluded by the lower bound on the Z' mass mz' = my 
resulting from the dedicated experimental searches |50 |. 


E. Halo stability 


If dark matter has significant self interactions and 
a suitable mass spectrum, the dark matter halos can 
collapse and reduce to dark disks via cooling through 
bremsstrahlung or Compton scattering on CMB photons. 
Due to the absence of a light dark photon, the feeble in¬ 
teractions and the mass scale of the MCP dark matter 
candidates we are proposing, we expect such a cooling 
process to be negligible in the present case. To show 
this, we assume the dark plasma is virialized and has a 
virial temperature 


^ _ Mmy 

-^vir — * \ ^2 


vir 


(37) 


where Udof = 3 is the number of degrees of freedom car¬ 
ried by a single particle and M and Rvir denote the mass 
and the radius of the virial cluster respectively. The char¬ 
acteristic timescale for dark bremsstrahlung cooling in 
our model is larger than the age of the Universe by many 
orders of magnitude m 


th 


rems 


3 

16 


'"‘V -‘^vir 

a^py e® 



(38) 


We proposed the most minimal renormalizable model 
of spin-1 MCP, based on the SM gauge group extension 
SU{2)l X U{1)y X SU(2)d, containing only a doublet 
Higgs field in the dark sector which is milli-charged under 
the U{1)y SM gauge group. We studied the basic phe¬ 
nomenological constraints of this model with the stress 
on a having a viable vectorial milli-charged (component) 
of dark matter. 

Constraints on the milli-charge ee versus the mass my 
were analysed by requiring that the vectorial MCP satis¬ 
fies the constraints imposed by the observed dark matter 
relic abundance. The production of the latter by both 
freeze-out and freeze-in has been investigated. In or¬ 
der to propose a first simplified scenario based on the 
proposed framework, we worked under a simplifying as¬ 
sumption that bars the effects of both the dark Higgs 
boson and the Z' component of SU(2)]j in determining 
the relic dark matter abundance. 

In this simplified scheme, we found that the freeze-out 
mechanism yields an overabundance of vector MCP with 
respect to the measured dark matter abundance. We be¬ 
lieve that relaxing our working assumptions could help to 
achieve the desired relic abundance within the freeze-out 
scenario and briefly commented on this possibility which 
we intend to investigate in a dedicated follow-up paper. 
On the other hand, within the freeze-in scenario, match¬ 
ing the dark matter relic abundance imposes e ^ 10 “^® 
regardless of the MCP spin. Owing to the peculiar be¬ 
haviour of the 77 —>■ VV cross section at high energy, 
spin -1 MCP present an additional production channel. 
In this case reproducing the dark matter abundance and 
respecting the mass constraints imposed by Z' searches 
however requires a reheating temperature of the order of 
Planck mass. 

For the values of mass and milli-charge emerging from the 
analysed freeze-in scenario, our MCP dark matter candi¬ 
date do not participate in the acoustic plasma oscillations 
at the recombination era and, therefore, avoid the severe 
bound cast by CMB analyses. In a similar fashion, the 
constraints emerging from dark matter halo stability are 
met owing to the smallness of the milli-charge e « 10 “^® 
required by the freeze-in mechanism. This value is com¬ 
parable with the condition for shockwaves formation in 










9 


cluster collisions that plasma physics indicates, possibly 
allowing for a test of the scenario owing to the implied 
astrophysical effects. 
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Appendix A: Cross sections and decay widths 

Here we report the relevant cross-sections at the lead¬ 
ing order in The cross-sections for // —)■ V~^V~ is 


Tff^V+V- 


24 


-jiyPf ^ det(s — M) ^ (l2my -I- 20TOyS -I- s^) 


(Al) 


X -f s(2m^Y^ - + Tr))- 

— s^(m^f{Yl — QY^Yfi + Y^) + 2m'^Yji{YL + Tr)) -1- s^(Y^ + Y^)^ , 


where (3x = a /1 — 4m^ /s and Yj^, Yji denote the hypercharge of the left- and right-handed fermion correspondingly. 
The determinant entering the propagators contains the mass matrix explicitly 


det(s - M^) = s((s - mlj(s - m|,) - - m^)) = s{s - m\){s - m\,). 

The cross-sections for W~^]V~ —> V~^V~ and 77 —> V~^V~ are 


2 2 

Tra € 


aiv+iv-^v+v- = 432^4 det(s - M) ^ x s {Urriy + 20myS + s^) 

X ^12774^ 4- 20tti^s , 

~Pv ^ 4 g ~ l^v) X (2 — Py — (1 — /?y)/3y^atanh(/3y)) -I- . 


77->y+v- 


In the following we assume that mw^ruf my and 
thereby neglect the SM masses. The high energy asymp¬ 
totic of the cross-sections reads 


(A2) 


(A3) 


(A4) 


(Tw+W-^V+V - 

(Jff^V+V- = 


432c((,s 

2 2 

Tra e 

432c^s 

STra^e^ 

TO?, 


X 18{Yl + Yi), 


and at the threshold s « 4my 


^w+w~->-v+v - 

aff^v+v- = 
^ 77 — 


Tra^e^ Py 
lUci 


2 > 


w 


(A5) 

(A 6 ) 

(A7) 

(A 8 ) 


r 


The asymptotics for the VV annihilations into the SM 
particles are 


327ra^e‘* i Qira^e^ 






9m y 


16ci 


S-l -bC>(s- 2 ). 


(All) 


Although the f-channel or term is of higher order in 
the milli-charge, it dominates the asymptotic regime with 
large s or small f3y. 

The decay rate of the Z' (assuming mz' ^ mf,mw) 
to SM particles is 


(A8) 




r !Z'- 


.VK+ry- 


^^amy 

48c2 

^^amy 

48c?„ 


X 8(r2 + Yl), 


(A12a) 

(A12b) 
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Assuming a heavy Higgs and 2mv > mz' these are also 
the only decay channels. The Z' width is then 

r. = (A13) 

Appendix B: The Feynman rules for the milli-charge 
expansion 

In the perturbative prescription, with the mixing as 
expansion parameter, the SM Feynman rules remain un¬ 
changed while there are additional 2-leg vertices from 


mass mixing connecting the dark and visible sector: 



z'-i, 


Z' - z. 


7-7, 


z -1, 


z-z. 


Besides a possible Higgs portal, the above interactions 
are the only ones that connect the two sectors at a dia¬ 
gram level. The dark sector Feynman rules are the usual 
rules for a spontaneously broken SU(2). 
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